(N 

o 

(N 



(N 



C/2 



o 

(N 

in 

m 

p 

rn 
O 
(N 



X 



Quantum phase transitions to superfluid state of chains in a polarized gas of 

dipolar molecules 
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We analyze the nature of quantum phase transition to a superfluid state of flexible chains in a gas 
of polar bosonic molecules confined in a stack of A'' identical Id ("cigar" type) optical lattice layers 
and polarized perpendicularly to the layers. Monte Carlo simulations within the J-current model 
show that, in the absence of the inter-tube tunneling, the transition to the A'^-layered superfluid is 
in the Berezinskii-Kosterlitz-Thouless universality class in the one-particle density matrix channel. 
The inter-layer tunneling changes it to the q = N 2d Potts universality. The low energy field 
descriptions of the transition are discussed in terms of conformal field theories. 

PACS numbers: 67.80. bd, 67.80.-s, 05.30.Jp, 61.72.Ff 



I. INTRODUCTION 

Quantum properties of extended objects ~ high 
energy strings (see in Ref.[l|), stripes in high-T^ 
superconductors |2| , vortices in superfluids and su- 
perconductors, dislocations in quantum crystals, 
etc. - arc of great interest to many areas of physics. 
Recent breakthroughs in creating and trapping 
high density samples of (polar) molecules [3] open 
up the possibility of realizing quantum chains of 
polar molecules. 

Self-assembly of classical chains has been inves- 
tigated numerically in Ref. [J]. So far, quantum 
chains have been studied in various approxima- 
tions which neglect tunneling of particles along the 
chains. In Ref.Q it has been proposed that stiff 
dipolar non-interacting quantum chains should 
form Bose-Einstein condensate. Inter-layer pair- 
ing in bilayred dipolar systems has been studied in 
Refs.[6|. Chains of fermionic molecules in multi- 
layered systems have been discussed in Refs.Q- 

Chains of indistinguishable particles as quan- 
tum objects represent a strongly interacting sys- 
tem which, in general, is not amenable to the stan- 
dard mean field or perturbation expansion meth- 
ods. Even the second quantization for chains is 
generically very problematic to implement. The 
complexity comes from the interplay between in- 
ternal degrees of freedom, that is, deformations of 
each chain and their centers of mass. Furthermore, 
several chains can exchange their segments with 
each other. Thus, on one side, identity of each 
chain is ill-defined, and, on the other, each chain 
looks slightly different from others due to the in- 
ternal excitations M- Under these circumstances 



it is not clear what type of order quantum chains 
can form. 

A numerical study of quantum chains which 
takes into account the partial-chain exchanges as 
well as the intra-chain dynamics has been per- 
formed in Ref. Q in the approximation of zero 
inter-layer tunneling. It has been shown that polar 
molecules in the iV-laycred geometry can form flex- 
ible (quantum rough) chains, and these chains can 
undergo a quantum phase transition to a superfluid 
phase characterized by off-diagonal long-range or- 
der (ODLRO) in the iV-body density matrix, while 
all M-body density matrices with M < N show 
insulating behavior (typical for Mott insulator) re- 
gardless of the filling factor provided it is the same 
in each layer. If the inter-layer (dipolar) inter- 
actions are weak, a stack of N layers features a 
TV-component superfluid (N-SF). Once the inter- 
action becomes stronger, the non-dissipative drag 
between layers will eventually convert the N-SF 
phase to the flexible chain superfluid (CSF) char- 
acterized by ODLRO only in the TV-body density 
matrix [9[. The corresponding transition is con- 
tinuous in the Id-geometry and discontinuous in 
the 2d-geometry for iV > 2 0]. The exact nature 
of the transition in Id was not, however, precisely 
demonstrated. 

In the present work we analyze how a finite inter- 
layer tunneling affects the transition. Since in the 
2d-layers case (iV > 2) the transition is of 1st order 
and weak tunneling cannot change this, we concen- 
trate on the Id-geometry (see Fig[T]), that is, when 
layers should rather be called as tubes. 

Our main findings arc as follows: The quantum 
phase transition into the iV-chain superfluid is in 
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FIG. 1: (Color online) Schematics of Id parallel to 
each other (three) layers each containing on average 
equal number (four) of dipole particles — circles with 
arrows depicting polarization perpendicular to the lay- 
ers (along Z-axis). Each layer represents Id optical 
lattice along X-axis. Tunneling along Z-direction is 
allowed, so that, while on average each layer contains 
the same number of particles, the relative numbers can 
fluctuate. The attractive interaction is mostly along 
Z-axis. 



the universality oi q = N 2d Potts model. That 
is, for TV = 2,3,4 the transition is continuous and 
for A'' > 4 it is of first order. For N = 2 the ex- 
plicit second quantized description is presented in 
terms of two-component Majorana fcrmions and it 
has been conjectured that the cases A^ = 3, 4 can 
be described by, respectively, 0(6) and 0(8) sym- 
metric conformal field theories. The numerically 
found i^-exponcnt (of the correlation length) val- 
ues are consistent with the analytical results for 
N = 2,3,4. It is confirmed that the transition 
becomes of 1st order for A^ > 4. It is also ex- 
plicitly shown that the transition to CSF at zero 
tunneling is in the Berezinskii-Kosterlitz-Thoulcss 
(BKT) universality class. 

In the next section, we will give an argument 
that the transition should be in the 2d Potts uni- 
versality. Then, in Sec. IIIII we will describe the 
microscopic Hamiltonian responsible for the chain 
formation and will introduce the coarse-grained 
dual version of the Hamiltonian in the discretized 
time approximation. The quantities used to char- 
acterize the transition and which are evaluated nu- 
merically will be introduced. The main numerical 
results for the correlation length criticality will be 
presented as well. Finally, the effective field model 
will be discussed in Sec. HVl 



II. ORDER PARAMETERS AND 

UNIVERSALITY OF THE TRANSITION 

TO THE CSF PHASE 

In the case of finite inter-layer tunneling (along 
Z in Fig. [T]), the one-particle superfluid (SF) is 



characterized by a single phase (p{x,z) which d(= 
fines the one-particle order parameter {'ip{x, z)) ~ 
e^'^, where 'ip{x,z) stands for the bosonic operator 
and in Id its mean is understood in terms of the 
algebraic order. In the CSF phase such order pa- 
rameter becomes zero and the condensate remains 
only for the product 
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ip{x,z). 
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So, the transition occurs between two supcrfluids: 
one characterized by (ip) ^ 0, (*) ^ (SF) and 
the other by (V') = 0, (*) ^ (CSF) (where the 
ODLRO in Id is also understood in the algebraic 
sense). In the next section we will give more ac- 
curate definition in terms of the A/-body density 
matrices. 

It is important that ^, Eq.(IT]), is invariant 
with respect to the transformation ij^^x, z) — > 
exp(27rzm(z, x)/N)ip(x, z) where to(z, x) is defined 
modulo N and obeys the constraint ^^ to(z, x) = 
pN,p = 0,1,2,.... Thus, m{z,x) can be broken 
as m = m' + rh into the discrete global part 
rh — p, p — 1,2,...,N — 1 and the local gauge- 
type m'{z,x) obeying X^z"^' = 0- Setting aside 
the discussion of a possible role of the local-gauge 
symmetry, we note that this discrete global trans- 
formation corresponds to the Potts model symme- 
try. Thus, it is natural to anticipate the q ~ N 
Potts universality for the quantum SF-CSF tran- 
sition. Accordingly, for Id tubes it should be con- 
tinuous (2d Potts) for N = 2,3,4 and discontin- 
uous for A^ > 4. The case N = 2 corresponds 
effectively to pairing transition of one-component 
bosons. Such transition has been proposed to be 
in the Ising universality (which coincides with the 
q = 2 Potts class) in Ref.[l3]. 



HI. HAMILTONIAN AND ITS 
J-CURRENT VERSION 

The microscopic Hamiltonian H describing SF 
and CSF is formulated in terms of the creation- 
annihilation operators aj.^,, axz of a boson at site 
X belonging to zth layer: 

^ / ^ ^z,z'^xz^^' z' '"^ / ^ * xz\x' z''^xz'^x' z' •/ 

{xx'),z,z' xz;x'z' 

(2) 
where tz^z' stands for a matrix of tunneling ampli- 
tudes: for z = z' ii is the intra-layer tunneling am- 
plitude til between nearest-neighbors, x' = x ±1, 



and for z' 7^ z it describes the inter-layer tun- 
neling i_L between sites x = x' in the neighbor- 
ing layers, z' = z ± 1 [no tunneling between dif- 
ferent sites located in different tubes is consid- 
ered]; Uxz = o.tz'^xz denotes onsite density oper- 
ator obeying the hard-core constraint; V^z-x'z' de- 
scribes the matrix clement for dipole-dipole inter- 
action between sites (xz) and (x'z'). It is charac- 
terized by strength Vd = dl/b^, where dz stands 
for the induced dipole moment and bz denotes dis- 
tance between two nearest layers. This interaction 
is mainly attractive along Z and repulsive along 
X. 



A. M-body density matrix 

The order can be characterized, in general, by 
the M-body density matrix 
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where 
aging. 



stands for the quantum-thermal aver- 



In the Id SF Di(x, z; a;', z') - l/lx- 



o.'l'' 



6< 1, 



exhibits algebraic order at large |a; — x'j. In the 
CSF, Di{x, z; x', z') ^ exp(-|x - x'|/Co), Co ^ 1, 
that is, it becomes short ranged at T = regard- 
less of the filling factor. Despite that, the A^-body 
matrix exhibits the algebraic order. On one hand, 
this order is characterized by the exponential decay 



Dn{xi 



,x-^, ..., Xjy^j '^ cxpl^ |a^mi 



J/Co) 



with respect to any pair of coordinates from ei- 
ther the set xi,...,Xm or x[,...,x[^. On the other 
hand, there is the algebraic order D^ ~ 1/| i?cm — 
^cml'^i c < 1, with respect to the "center of mass" 
coordinates Rem = [xi + ... + xn]/N and R[^ = 
[x[ + ... + x'j^]/N defined, respectively, for the first 
and the second sets of the coordinates, provided 
\Rcm - Xm.\ < ■Co and \R'em - x'^\ < Co for aU to. 

It is important to note that, while all Dm, ^I = 
1,...,N - 1 exhibit algebraic ODLRO in the SF 
phase and exponential decay in the CSF, D^ is 
long-ranged in both phases — SF and CSF — with 
respect to \Rcm — R'cml ~^ '^ (^^ the algebraic 
sense). Thus, the transition SF to CSF can be de- 
tected by critical behavior of any Dm, M < N. 
We also note that the same criticality controls 
the long-distance behavior of Dn with respect to 
I -Rem - Xjn\ (or \Rem - x'^\). That is, in SF 
phase Dpj is trivially long-ranged with respect to 
I -Rem — Xm\ because ZJjv can simply be factorized 
into a product of Di. In contrast, in the CSF- 
phase, while exhibiting ODLRO with respect to 



00 


-cKGjWdU 




^/ 


80 


- L=300 




0.60 0.61 0.62 0.63 U 



0.600 0.605 0.610 0.615 0.620 0.625 0.630 0.635 U 

FIG. 2: (Color online) d{Gz) /dU versus the interaction 
strength f/ for L = 60, 100, 200, 300 with P = L. Inset: 
(Ga) versus U for L — 100. The SF phase corresponds 
to (Gs) fa 1 and the CSF to (G3) « 0. The transi- 
tion point SF-CSF for a given size can be identified by 
the maximum of dG^/dU reaching the thermodynam- 
ics limit at Uc ~ 0.61. 



Rem ~ -Rem' ^N IS short-rangcd with respect to 
I -Rem - Xm\ (or |i?^„ - x',J) . Thus, the criticality 
can also be detected by measuring the behavior 
of the relative distances Xm (or x'^). Specifically, 
we have considered the square of so called gyration 
radius [9| as the mean of 



R^g - 



- y 

m,n— 1,2,. 



(4) 



,N 



with respect to the first set of the coordinates of 
Dn, provided the coordinates from the second set 
are kept within some distance ~ Co from R'^^ 11 il l . 



In the SF of a length L, R^g ^ 
^i^L^ and in the CSF i?|p^ 



Rl ^ 0{L^) « 
^ 0(1) « Co. In 
what follows we will be calculating the mean of 
the ratio Gn = R^/Rq, so that it is changing from 
Gat « 1 in the SF state to Gat ~ l/L^ r^ in 
the CSF phase. It is worth mentioning that Rg 
can be viewed as a typical width of a chain. For 
strongly bound case this width is ~ Co ~ li and 
in the SF phase it is ^ L, and, thus, it exhibits 
critical behavior typical for correlation length. 



B. J-current formulation 

Hamiltonian ([2|) and the gyration radius (|4]) have 
been used for ah initio simulations of a single chain 
(with exactly one polar particle per layer) for the 
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FIG. 3: (Color online) d^Gz) /dU versus (G3) for sizes 
L = 20,40, ..., 300 rescaled by a factor \{L) in order 
to achieve collapse to the curve L = 100 (A(IOO) = 1). 
Inset: d{G-i) /dU versus (G3) for the same sizes. 



FIG. 6: (Color online) The maximum value of 
d{GxY)/dU versus InL for the case of zero inter- 
layer tunneling. A*' = 4, K{z) = 00. The solid red 
line is the fit by the finite size scaling ansatz for 
the BKT-transition: d{GxY)/dU = Aln^(L/Lo), A = 
0.205, Lo = 1.98. 
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FIG. 4: (Color online) The rescaling factor \~^(V) ver- 
sus L for N = 3 from Fig [3] The slope gives the corre- 
lation length exponent v — 0.835 ± 0.015. 
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FIG. 5: (Color online) Bimodal histogram of energy 
P{E) for iV = 5 tubes with L = 400, 13 = 400. The 
first-order transition happens at Uc = 0.7235. 



case i_L = 1^]. It was found that the chain can 
undergo quantum roughening transition with the 
tuning parameter being the interaction strength 
Vd- The transition is, practically, insensitive to 
the interaction range. 

The simulations at finite densities n in each 
layer have been conducted in the discrete-time J- 
current-type formulation [12[ of the Hamiltonian 
([U Q- For the purpose of analyzing the uni- 
versality of the transition this approach turns out 
to be much more efficient than the ab initio one. 
Here we will be using similar model where the 
inter-layer tunneling is allowed. The actual dipole- 
dipole interaction will be replaced by onsite attrac- 
tion between neighboring layers, with the intra- 
layer dipole-dipole repulsion ignored. The corre- 
sponding space-time action, then, becomes 



^^ = E 



K{b){Jt)^ , U{b)iV,J,y 



-M 



(^') 



(5) 

where Jb is the integer bond current obeying Kirch- 
hoff's conservation law [12]; the summation is per- 
formed over all space-time bonds b (coming out 
from a space-time site (x, r, z) either along ±x or 
along imaginary time ±f or along ±z directions); 
VjJfe = Jii{x,T, z + l) — Jii{x,T,z); ^ denotes chem- 
ical potential. [Here we tuned /i to have 1/2 filling 
of bosons per site in each tube] . Periodic boundary 
conditions along space < a; < £—1, < z < iV— 1 
and along imaginary time < r < /5, with 13 ~ L, 



where L ~ 2,3,...., have been used. The coef- 
ficients K, U can be related to tz,z' and Vxzyx'z' 
from Eq.©: K{z) ^ l/t±, K{x) = K{t) « l/i||, 
U ~ Vd- The case K{z) ~ oo corresponds to zero 
intcr-layer tunneling (studied in Rcf. [0] ) . Here we 
will focus on K{z) = K{x) = K{f) situation as 
the one which naturally represents the whole uni- 
versality class. 



We note that the action ([S]) can be viewed as a 
coarse grained dual representation of the Hamilto- 
nian ^. While being not precise for quantifying 
finite energy (non-universal) properties of the sys- 
tem, the J-current model [12| belongs to the same 
universality class as the original model ^. Thus, 
for the purpose of this work and for sake of numer- 
ical practicality, it will be sufficient to study the 
model ([5]). 



We have performed Monte-Carlo simulations 
of the model ([S]) within the Worm Algorithm 
approach [l3|. Green's function in imaginary time 
(as well as the density matrix ([3|)) is given by 
the statistics DAr({xm,rm, ^m}; {a;™, t^, ^m}) of 
"sources" and "sinks" of the bond currents located, 
respectively, at {xm,Tm,Zm), rn = l,2,...,iV, and 
(x^, r^, z!^), m= 1,2, ...,N, lattice points. In or- 
der to insure the condition |i?cm ~ ^'m\ — ^o, 
while {xrn,Tm,Zrn), m = l,2,...,N, are 
free to take any value, we have convo- 
luted DNi{Xm,Tm,Zjn};{x'^,T^,z'^J) with 

P = CXp{~j:rn,n[\^'m ~ <\ + K. ^ 

Tn\]/£.o) as DN{{Xm,T^,Zm};Rcni) = 

J Dx'Dt'Dz'DnPS [R',^ - E„ x'JN) 
and, accordingly have evaluated the means 
of the normalized gyration radius (Gat) and 
of the center of mass distance {\Rcm — R'cm\) 
where (...) = Z-^ J OxDrDzdR'^^.-DN, Z = 
J DxDTDzdR'„-,^DN. 



For sake of numerical efficiency we have sym- 
metrized the model (O by choosing U{b) indepen- 
dent of the type of a bond, that is, U{b) = U. 
The CSF phase has been identified by the condi- 
tion (|i?cm ~ R'cm\)/L = const and {Gn) ^ o(L~^) 
for U > Uc, where Uc corresponds to the quan- 
tum critical point (QCP). In the SF phase (that is, 
U < Uc), while the first condition remained, prac- 
tically, unchanged, (Gat) « 1 with high accuracy. 
The criticality of the SF-CSF transition has been 
analyzed through evaluating the divergent behav- 
ior of d{GN)/dU in the vicinity of C/ = Uc- 



C. Finite size scaling of the gyration radius 

As discussed above, d{Gt^)/dU exhibits singu- 
larity in the limit L — >■ oo. The change from 
(Gjv) w 1 to {Gn) ~ occurs in a narrow range 
SU = \U — Uc\ around the critical point Uc ^ I- 
Such a behavior is clearly seen in Fig. [2J the range 
of the transition 5U narrows as L increases. 

This range is controlled by the diverging cor- 
relation length ^{U) - |C/ - Ucl^'', ly > 0, 
where ly stands for the correlation length expo- 
nent. According to the finite size scaling approach, 
(Gjv) can be represented as some regular function 
P{y)^ V = L/S,{U) varying from F{y = 0) = 1 
to F{y = oo) = over the range y ~ 1. Thus, 
d{GN)/dU w F'y/6U - L^/^ Loosely speak- 
ing, one can view this relation as d{GN)/dU ~ 
1/SU, SU sa L-i/'' ^ 0. 

We have evaluated this derivative numerically 
by Monte Carlo [IJl and constructed the graphs 
d(GN)/dU versus (Gat) by scanning over U around 
the critical point Uc for sizes L = /3 = 10, 20, ...300. 
These graphs turn out to be self-similar so that 
d{GN)/dU for all sizes > 10 collapsed on a sin- 
gle master curve by simple rescaling of d{GN)/dU 
for size Li to another size L2 as d{Giq)/dU -^ 
\{L)d{G m) / dU . Then, the rescaling coeflicient 
A, which represents the inverse width 5U as A oc 
l/(5t/, have been plotted in the log-logL axes in 
order to determine the critical exponent v. The 
results of this procedure are presented on Figs. I3I4I 
for the case A^ = 3. [The same procedure has 
been used in the cases N = 2,4 as well]. The 
found exponents are: v = 0.972 ± 0.02 for N = 2, 

V = 0.835 ± 0.015 for A^ = 3, I' = 0.735 ± 0.015 
for A^ = 4. [The shown errors include statistical 
errors as well as the systematic errors due to the 
subleading contributions] . Wc note that the value 
oi V ioY N = 2 is consistent with the d = 2 Ising 
(or q = 2 Potts) universality. Wc also note that 
the values of v for A^ = 3 and A^ = 4 are consis- 
tent with the corresponding ones v = 0.837 and 

V = 0.756 obtained by the Renormalization Group 
calculations for the q = 3,4 2d Potts model |14| . 

For Af > 4, the transition was found to be of 
first order. It has been detected by the bimodality 
of energy histogram. Fig. [S] While for A^ = 5 such 
bimodality develops on sizes L > 400, /? = i, for 
A^ = 8 it is already well developed ai L — IQQ, ji — 
L. 

The finite size analysis has been applied to the 
case of zero inter-layer tunneling as well, when the 
transition is expected to be in the BKT univer- 
sality. That is, ^ ~ exp(...|C/ - Uc\-^''^). The 



variation of the gyration radius (Gxy) in. this 
case ean also be represented by some regular func- 
tion F(i/) characterized by the range y ^ 1 with 
y = L/i . Thus, d{GxY)/dU « F'y ^ [SUy^/^ ^ 
{\x\{L/Lo)Y (where Lq stands for some microscopic 
scale) at its maximum. The maximum value of 
this derivative has been plotted as a function of 
L = 30, ..., 600 in FiglHl As can be seen the fit of 
the data is consistent with the In L dependence 
with high accuracy. 



IV. THE EFFECTIVE MODEL IN TERMS 
OF MAJORANA FERMIONS FOR TV = 2 

Here we present the effective model describing 
the SF to CSF transition in the approximation 
where the dipole-dipole interaction is reduced to 
the attraction between nearest neighbors layers. 
Taking into account the single occupancy con- 
straint we can rewrite Hamiltonian ^ in terms 
of Pauli matrix operators: 



H - ^iN = 



E 

3 



{ E [-MMtz''3+^,z + h.C.) + Va(Tl^(7^^^^^^ 
z = l,2 

] + ti- {^s^Ja + h.c.) ~ Vial^al^ } (6) 



J,z 



where cr* operators stand for the bosonic cre- 
ation and annihilation ones a\a in Eq.Q and 
h = <T'- + 1. While in the SF phase both fields 
a^-^ and a~2 condense (in the algebraic sense), the 
CSF phase {N = 2) corresponds to the condensa- 



tion of ^(j) = o'7iO'^2' with a~^ being disorderecR 
Another way to say this is in terms of the phases 
(fiz of the fields a~^ ^ exp(i(p^). The backscatter- 
ing events can make ipi — ip2 strongly fluctuating 
while ifi + ip2 remains well defined. Accordingly, 
cos((y9i — ip2) becomes irrelevant |15| . This consti- 
tutes the transition point from SF to TV = 2 CSF. 

In fact, similar mechanism works in mixtures 
of two non-convertible atomic species [la] (where 
i_L = 0). It is important to note, however, that, 
despite such a similarity, the universality of the 
transition dramatically depends on t±: while for 
ij^ = it is of BKT-type [ig] , it becomes of Ising- 
type at t_L 7^0 ^. 

We will treat this model using bosonization tech- 
nique (see in Ref. [1^ ) which in the context of the 
model © was pioneered by Schulz ^^. For com- 
pleteness we will reproduce the calculations here. 

The Hamiltonian ^ for two independent tubes 
in the low-energy limit (ignoring the backscattcr- 
ing events) is the same as for the spin S=l/2 XXZ 
model. In the continuous limit it is equivalent to 
the Gaussian model 



Ha=^^ 



dx 



'K-\d,'^a?+K{d.,ea?], (7) 



where a = 1, 2 labels the tubes and 8a is the field 
dual to $„: [dxQaix),<^a{y)] = -iTTS{x - y). The 
Luttinger parameter K is determined by the intra- 
chain interactions. If for convenience we assume 
that til > then the continuum limit of the opera- 
tors is given by the following bosonization formu- 
lae: 



J 



<yt{x) 



„lV2ir*„ 



G 



,l\/27r(e„+*a)-K21fcFa; , lV27r(-e„+$„)- 



■2\kFX 



(27r)i/2 
al{x) = -^dxOa + TTTTTT^ smC^kpx + V2^Qa){-l) 



(8) 



where dots stand for less relevant operators and 
C, C^ are amplitudes determined by the short 
range physics. The Fermi momentum kp for each 
chain is determined by its chemical potential [i so 
that A:p(yU = 0) = 0. In the following we will always 
assume that the chemical potential is non-zero so 
that the spin fluctuations are incommensurate with 
the lattice. 



Substituting ^ into © and defining the fields 
$1,2 = \/^(ify'$+ ± i^y $_) , 



ei,2 = ^(Av^/^e+±A'i^/^e_ 

\dx^a{x),^h{y)'\ = -i(5a6(5(x - y). 



(9) 



we obtain the Hamiltonian H — H^ + H- , where 



H+ describes the symmetric mode (+): 



H, 



Ax 



(10) 



and H- contains only the anti-symmetric fields: 

F_ = ydx{^[(a,$_)2 + (a,e_)2]+ (ii) 

Vjcos(^AtiK^^\ - Vc cos f v/47r/A'_e_U, 
where Vc ^ Vi,Vj ^ t± and 



K±= K± 



2'KV 



(12) 



At this juncture we note that the Hamiltonian 
(fTO|) describes a mode which is not critical at the 
SF to CSF transition. In other words, it is the N = 
2 CSF. In contrast, the Hamiltonian (TTll) accounts 



J 



for the transition so that $_ 
in the CSF. 



becomes disordered 



Depending on which of the cosines in pT|) takes 
over, the ground state of this model describes 
either quasi long range superfluid order or pair 
density wave. The latter state has a singular- 
ity in the density-density correlation function at 
the finite wave vector 2kp. When both cosines 
are relevant (that is at 1/2 < K^ < 2) these 
states are separated by a QCP, the location of 
which is approximately determined by the relation 
(V>/A)^- - {Vc/Ay/^-, where A is the ultravi- 
olet cut-off. Fulfillment of the above condition on 
K- is essential for the subsequent arguments. 

The vicinity of the QCP can be studied analyt- 
ically when K^ « 1 (this will be our assumption 
throughout the rest of the paper) . In that case it 
is convenient to refermionize (fTTj) with the result 



{IV— 
-^{-PRdxPR + PhdxPL - ilBdxim + riLdxTlL) + 

4:TTV-{K-. - 1)prPlVb.i1l + 2im+pB.PL + 2i'm-riRriL>, 



(13) 



Where ■m± = Vj zLVc and pl,r and riL^n are left- 
and right-moving components of Majorana (real) 
fermions. This model is equivalent to the con- 
tinuum limit of two quantum Ising (QI) models 
coupled by the energy density operators [l^- The 
Monte Carlo resuh i' = 0.972 ± 0.02 for the model 
^, N ^ 2, which corresponds to the exact solu- 
tion i^ ~ 1 of d =: 2 Ising model, is consistent with 
such conclusion. 

In the cases N = 3, 4 wc conjecture the tran- 
sition to be described by, respectively, 0(6) and 
0(8) parafermion model perturbed by the energy 
density operator. The scaling dimension d of this 
operator is the same as in SUAr(2) Wess-Zumino- 
Novikov-Witten model: d = 4:/(N+2) and, as con- 
sequence, ly = 1/(2 - d)[23|,|2li. So, for iV = 3 we 
get d = 4/5 and i^ = 5/6 w 0.833 and for A^ = 4 
we get d = 2/3, that is, v = 0.75. These values 
are consistent with the above Monte Carlo results 



V = 0.735±0.015 for A^ = 4, and t/ = 0.835±0.015 
for A^ = 3. 

When this work was prepared for publication we 
learned about the preprint by Lecheminant and 
Nonne [22| which results have a substantial overlap 
with ours. 
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